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Abstract 



In PP a two-Higgs-doublet model with maximal symmetry under generalised CP trans- 
formations, the MCPM, has been proposed. The theory features a unique fermion mass 
spectrum which, although not describing nature precisely, provides a good approxima- 
tion. It also predicts the existence of five Higgs bosons with a particular signature. 

In this thesis I implemented the MCPM into the Monte-Carlo event generation package 
MadGraph, allowing the simulation of any MCPM tree-level process. The generated 
events are in a standardised format and can be used for further analysis with tools 
such as PYTHIA or GEANT, eventually leading to the comparison with experimental 
data and the exclusion or discovery of the theory. The implementation was successfully 
validated in different ways. It was then used for a first comparison of the MCPM signal 
events with the SM background and previous searches for new physics, hinting that the 
data expected at the LHC in the next years might provide exclusion limits or show 
signatures of this model. 



Zusammenfassung 

In m wurde ein Zwei-Higgs-Doublet-Modell mit maximaler Symmetrie unter generalisier- 
ten CP-Transformationen vorgeschlagen. Dieses sogenannte MCPM hat ein einzigartiges 
Massenspektrum, das die Natur zwar nicht exakt beschreibt, aber eine gute Naherung 
darstellt. Das Modell sagt die Existenz von fiinf Higgs-Bosonen mit speziellen Eigen- 
schaften voraus. 

In dieser Arbeit habe ich das MCPM in die Monte-Carlo-Simulation MadGraph imple- 
mentiert. Damit lassen sich Ereignisse fur beliebige Tree-Level-Prozesse generieren, die 
dank standardisierter Schnittstellen mit vielen Programmen wie PYTHIA oder GEANT 
weiter verarbeitet und schliefilich zum Ausschluss oder zur Entdeckung des Modells in 
Beschleuniger-Experimenten verwendet werden konnen. Die Implementierung wurde mit 
verschiedenen Methoden erfolgreich validiert. Anschliefiend wurde damit ein erster Ver- 
gleich einiger MCPM-Prozesse mit dem Untergrund aus Standardmodell-Prozessen sowie 
mit bisherigen experimentellen Studien durchgefiihrt. Die Ergebnisse deuten darauf hin, 
dass sich aus den Daten, die am LHC in den nachsten Jahren erwartet werden, entweder 
Ausschlussgrenzen oder Spuren des MCPM gewinnen lassen konnten. 
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1. Introduction 



The Standard Model of particle physics (SM) is one of the success stories of theoretical 
physics. This quantum field theory, which is based on a SU (3)c x SU (2)lxU (l)y gauge 
group, describes all the fundamental particles that have been observed so far and the 
strong, electromagnetic and weak interactions extremely well. Only the Higgs boson it 
predicts is yet to be discovered. 

However, from a theoretical point of view the SM is somewhat unsatisfactory. Its particle 
and mass structure seems arbitrary. There are three generations of fermions, but no 
explanation for this number. Of these fermions, the t, b and r are much heavier than 
the corresponding particles of the first and second generation^ 
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(1.1) 



There is no explanation for this mass hierarchy, either. Finally, the SM features the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix V. The absolute values of its entries are 
given by [2] 

f\V n \ | Via | |Vi 3 |\ /0.97 0.23 0.00\ 

|V 2 i| IV22I jVaa] « 0.23 0.97 0.04 , (1.2) 
\|V 3 i| |V 32 | \V33\J \0.01 0.04 1.00/ 

which is to a good approximation unity. Again, the SM does not provide any reason for 
this. Not counting the representation assignments, there are 18 free parameters in the 
theory [3]. 

A lot of research has been done on possible extensions of the SM. One ansatz is the 
extension of its scalar sector from one scalar doublet to two scalar doublets. By requiring 
that such a two-Higgs-doublet model (THDM) is invariant under so-called generalised 
CP transformations, a model with intriguing features can be constructed: the maximally 

1 Here the MS quark masses evaluated at fj, = vq ~ 246 GeV have been used [see[T]. 
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1. Introduction 



CP-symmetric two-Higgs-doublet model (MCPM). It has been developed and studied by 
various authors in the last years [H HHTU] . In the MCPM, the fermions of the first two 
generations are massless, while the t, b and r get a mass through electroweak symmetry 
breaking p]. The CKM matrix is the unity matrix. Although this is not precisely as 
observed, it is a good approximation to nature. 

The MCPM predicts the existence of five Higgs bosons. These couple to the fermions in 
a unique way. One Higgs boson behaves similar to the Higgs boson in the SM. However, 
four of the five Higgs particles couple only to the second generation fermions, although 
these are massless in the theory [TJ [5]. Thus the model has a phenomenology that is 
different from the SM and other extensions of the SM. In chapter [2] the MCPM and its 
phenomology are developed and explained. 

One obvious question remains unanswered: can this theory be either discovered or ex- 
cluded at the Large Hadron Collider (LHC) at CERN? To answer this, the experimental 
data has to be confronted to Monte-Carlo simulation results. In this Bachelor thesis, 
I implemented such a simulation of the MCPM based on the event generation pack- 
age MadGraph [11 . With this tool it is now easy to numerically calculate tree-level 
cross-sections for arbitrary processes and to generate events for further analysis. This 
constitutes the main part of this project and is documented in chapter |3j 

The MadGraph implementation of the MCPM was systematically validated in different 
ways. The results are presented in chapter [4j Using this new tool, the most relevant 
MCPM signatures were compared to the background from SM processes. It was also 
checked whether the results of the search for a Z' boson are likely to have already 
excluded the MCPM or a part of its parameter space. This analysis is the topic of 
chapter [5] 
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2. The maximally CP-symmetric 

two-Higgs-doublet model (MCPM) 



In this chapter the maximally CP-symmetric two-Higgs-doublet model is developed. At 
first, the framework of two-Higgs doublet models is introduced. Generalised CP trans- 
formations are defined and it is analysed under which conditions a THDM is invariant 
under these transformations. Based on these criteria, the MCPM is constructed. Its par- 
ticle spectrum and properties after electroweak symmetry breaking are derived. Finally, 
the phenomenology of the MCPM particles at colliders such as the LHC is sketched, 
including Higgs boson production mechanisms, their decay modes and the background 
from SM processes. 

2.1. General two-Higgs-doublet model 

A general two-Higgs-doublet model is a theory based closely on the SM. It features the 
gauge group SU(3)c x SU(2)l x U(1)y with gauge fields and B^, corresponding 

to colour, weak isospin and weak hypercharge, respectively. There are also the usual 
fermions: for each generation m there are up- and down-type quarks u m ,d m as well as 
charged leptons e m and neutrinos v m . The theory is chiral, i.e. left- and right-handed 
fermion fields are treated differently. The left-handed fermions are organised in SU(2)l 
doublets 



while the right-handed fermions itjj , d 1 ^, transform as singlets. Let us for now consider 
n generations of these fermions (later the observed three families will be assumed) . 

The gauge and fermionic sectors of the Lagrangian are given by 




and 




(2.1) 



^Fermions = qfilpq? + + + 4^ + 



R 



(2.2) 
(2.3) 
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2. The maximally CP-symmetric two-Higgs-doublet model (MCPM) 



following the usual conventions. A detailed explanation of the symbols and quantities 
used here can be found in appendix |A.l 



All of this is just as in the SM. However, the scalar sector is different. Instead of one 
complex scalar doublet, there are two such fields, denoted by 

^ (x)= (^(5) (z=i ' 2) - (2 - 4) 



Both of these scalar doublets are assigned weak hyperchargdj y = \ . Their properties 
are given by the scalar sector of the Lagrangian density, 

^Scalar = {D^(D^i) - V(ip U ip 2 ). (2.5) 

V(^pi,(p2) is the Higgs potential. Its most general gauge-invariant and renormalisable 
form is 

V(<pi,(p 2 ) = aijflfj + hlmn{<Pk ( Pl)( l Ptn l Pn) (2.6) 

with arbitrary coefficients aij,bki mn S 

The last part of the Lagrangian describing the THDM are the Yukawa terms coupling 
the fermionic to the scalar fields: 



-i- Yukawa — ~ 6rOj t Pi L L + u R L 'u ViQL a 0R L 'd { PiQ.L + tl - c - \ Z -<) 

The Yukawa couplings C* , C^, C\ are matrices in the generation space and 



~t T 



V 
-1 0, 



Now all parts of the THDM can be put together. 

Definition 1 (THDM) A general two-Higgs-doublet model is a quantum field theory 
described by the Lagrangian density 

&THDM = C-Gauge + £-Fermions + ^-Scalar + ^-Yukawa (2-8) 



where the different sectors are given in (2.3), (2.5) and (2.7). 



1 In the context of supersymmetric models, the two scalar doublets are usually assigned different weak 
hypercharges of | and — § respectively. However, this is only a matter of convention, see [I] for 
details. 
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2.2. Developing the MCPM 

2.2.1. Generalised CP transformations 

In quantum field theories charge conjugation (C) and parity transformation (P) play 
important roles. Especially interesting is their combination. The CP transformation was 
long thought to be a fundamental symmetry of nature, until its violation was discovered 
in 1964 [T2]. It is strongly related to the question why there is more matter than 
antimatter in the universe, which has still not been answered satisfactorily |13j . 

Under a standard CP transformation (denoted by CP S ) scalar fields transform as 

CP 

ifi{x) % tiixT) (2.9) 

and spinor fields as 
C P 

iP m (x) ^7° 7 y m V) (2-10) 



with 




°\ CPs , fx°\ 

>J= J (2-11) 



-x , 



due to the parity transformation |14[ section 3. 60 

In this standard CP transformation the spinor fields of the different generations do not 
mix, neither are the two Higgs fields transformed into each other. However, it is possible 
to define generalised CP transformations that allow these kinds of mixings [15 . For 
instance, such a transformation might consist of the usual charge conjugation and parity 
transformation, but additionally transform ip\ into if 2 and if 2 into ip±. 

In other words, generalised CP transformation are standard CP transformations together 
with a change of basis in the scalar doublet space and a change of basis in the generation 
space. These basis changes are not totally arbitrary. Just like two successive standard 
CP transformations never change the physics of a system, two successive generalised CP 
transformations are not allowed to change the physical states of the particles either. 

Definition 2 (Generalised CP transformation) A generalised CP transformation 
CP g transforms the spatial coordinates as in ( 2.11[ ), scalar fields tpi (i = 1,2) as 



CP g 



ipi(x) > U 9ij if*Ax') (2.12) 



2 The fields may gain additional minus signs under a CP S transformation. This is a matter of convention. 
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2. The maximally CP-symmetric two-Higgs-doublet model (MCPM) 



and spinor fields ip m with n generations as 
c P 

4> m (x) %iU^°j 2 i> k *(x'). (2.13) 

The matrices U v £ U{2), £ U(n) have to satisfy the condition that two successive 
gauge transformations do not change the fields at all or are at most equal to a gauge 
transformation: 

CPgOCPg^l. (2.14) 

The transformation properties of the gauge fields are described in IJjj. 

Generalised CP transformations can be classified and further analysed, see for instance 

HI- 

2.2.2. Requirements for the MCPM 

It is now interesting to analyse which conditions a THDM has to satisfy to be invariant 
under CP g transformations. In [1J the concept of maximal CP symmetry has been 
introducec 

Definition 3 (Maximal CP symmetry) A THDM is called maximally CP-symme- 
tric if for each U v G {1, o~ , ia 2 , <r 3 } there is a suitable S U(n) for each fermionic 
field such that the theory is invariant under the generalised CP transformation described 
by Utp and the different U^. Here a 1 are the Pauli matrices. 

It has been found [T] that this requirement is not enough to restrict the THDM to a 
theory with unambiguous mass spectrum and properties. Instead, it yields different 
possible scenarios. However, some of these cases predict a phenomenology that has 
already been ruled out by observations and can be excluded for this reason. 

A THDM describing nature has to fulfill the following conditions [T]: 

• When the scalar fields acquire vacuum expectation values, the electroweak 
SU(2)l x U(l)y symmetry is spontaneously broken to the electromagnetic U(1)q 
symmetry. The Higgs potential is stable in the sense of |6j. 



3 Despite its name, maximal CP symmetry does not require a symmetry under all generalised CP 
transformations. Requiring such a broad symmetry would necessarily lead to massless Higgs bosons. 
See [T] for more details. 

4 Note that in this thesis the abbreviation MCPM (for maximally CP-symmetric two-Higgs-doublet 
model) will exclusively refer to a specific model which will be defined later, not to any THDM with 
maximal CP symmetry. 
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• After symmetry breaking, the theory predicts none of the following: 

— Large lepton flavour-changing neutral currents. 

— Mass-degenerate massive fermions. 

— Massless Higgs bosons. 



2.2.3. Consequences 

The requirement of maximal CP symmetry together with the phenomenological restric- 



tions given above confines the Higgs potential (2.6) to a very specific form jU 



Theorem 1 The Higgs potential of a maximally CP-symmetric THDM which satisfies 
the criteria given in section \2.2.2\ has the form 

+ /j, 2 (-i<p\ip2 + ivl<Pi) 2 + M^i^i ~ ^2^2) 2 (2.15) 
with real parameters satisfying 



Ml > (Jq > 1^3, 


(2.16) 


?7oo > 0, 


(2.17) 


/"a + ??oo>0 for a = 1,2, 3, 


(2.18) 


£0 <o, 


(2.19) 


M3 < 0. 


(2.20) 



The consequences for the Yukawa couplings depend on the number of generations. The 
result for n = 1 family of fermions is surprising [1 : 

Theorem 2 In a maximally CP-symmetric THDM with one generation that satisfies 



the criteria given in section 2.2.2, all Yukawa couplings C^ nn , CJ 71 ™, Cf mn in (2.7) 
must be zero. 

In other words, in the maximally CP-symmetric THDM with only one fermion family 
no fermion will get a mass through electroweak symmetry breaking. This is remarkable: 
While in the SM the number of generations is totally arbitrary, the MCPM requires at 
least two generations for fermion masses to appear. 

For n = 2 generations the situation is different [T]: 
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2. The maximally CP-symmetric two-Higgs-doublet model (MCPM) 



Theorem 3 The Yukawa sector of a maximally CP-symmetric THDM with two gen- 
erations (numbered 2 and 3 for reasons that will become clear later) which satisfies the 



criteria given in section 2.2.2 is given by 



■ Yukawa 



-q (e%<p\l 3 L - 44^1) 

(— S T • 2 3 — 2 T. 9 2\ 
urVi W Ql ~ u R <p 2 w Ql ) 

- c d (d^Ul - d&tti) + h. c. (2.21) 



with real parameters c\,c u , Cd > 0. 



Here <p% couples to one generation of fermions (labelled 3), (f2 to the other (labelled 
2), both with the same coupling constants. This will lead to a specific mass spectrum 
after symmetry breaking: one generation will become massive, while the other remains 
massless, as will be demonstrated in section |2.2.4 So for two fermion families, the 



requirement of maximal CP symmetry yields an intrinsic mass hierarchy. 

The general case of n = 3 generations has not been fully analysed so far. However, a 
maximally CP-symmetric THDM with three generations, the so-called MCPM, can be 
created in a straightforward way: By identifying the solution for one generation with 
the first family (u, d, e, v e ) and the solution for two generations with the second (c, s, 
fi~ , Up) and third (t, b, r, u T ) families^] and combining them, the MCPM is defined [lj. 
In other words, the three-generation MCPM is introduced with the following properties: 
The second and third generation of fermions couple to the second and first Higgs doublet 
respectively, and the first family of fermions is uncoupled to the scalar fields. 



Definition 4 (MCPM) The THDM with three fermion generations, Higgs potential 



(2.15) and Yukawa sector 



(3 1" 3 2 "f - 2 i 
(— 1 T • 2 3 -2 T ■ 2 2\ 

+ Cu {URiPi la q L - u R (p 2 io q L J 

- c d (dy\ql - dy 2 q 2 L ) + h. c. (2.22) 



with parameters /X2, ^3) ^00) £o> Q, c u , c d E M satisfying (2.16) to (2.20) and ci,c u ,Cd > 
is called MCPM. It is maximally CP-symmetric, as can be seen from theorem and 
theorem 



3 Now the labels 2 and 3 for the two generations in ( |2.21[ ) make sense. 
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2.2. Developing the MCPM 



This MCPM is the topic of this thesis. I should stress that it does not follow necessarily 
from the requirement of invariance under generalised CP transformations. It is not the 
most general maximally CP-symmetric THDM with three generations. However, it is 
a very manifest choice and leads to interesting properties after spontaneous symmetry 
breaking, as will be seen soon. 



2.2.4. Electroweak symmetry breaking 



The Higgs potential (2.15) is minimised by vacuum expectation values (vev) for the 
scalar doublets. Without loss of generality these values can be chosen such that only 
the second component of (pi gets a non-zero vev v$ [T]: 

M = ^ (°) • (2.23) 
(**> = (fj , (2.24) 

where 



v = J (2.25) 

V »700 + A*3 

plays the same role as the vev in the SM. 

The scalar doublets can then be expanded around this vacuum state. In unitary gauge 
this reads [6] 

Mx) = 4= ( ,°// J , (2-26) 
V2 \vo + p (x) y 

^ = \^{h'{x) + \h"{x))) ■ ( 2 - 27 ) 
Here p'(x), h'(x) and h"(x) are real scalar fields, H + {x) is a complex scalar field with 

H~{x) := (h + (x)Y . (2.28) 

These fields are excitations from the vacuum state, so they represent physical particles: 
the physical Higgs bosons of the MCPM. The p' , h! and h" are uncharged and their own 
antiparticles, the pair turns out to be charged. As scalars all Higgs bosons have spin 
0. 



17 



2. The maximally CP-symmetric two-Higgs-doublet model (MCPM) 



The gauge fields behave exactly as in the SM. The electroweak SU(2) l X U{V)y symmetry 
is broken down to the electromagnetic U(1)q symmetry. The and fields form the 
massive and bosons and the massless photon [see for instance HH chapter 
20.2]. 



In addition, the scalar sector (2.5) with Higgs potential (2.15) yields mass terms for the 



Higgs bosons. These masses are given by jT] 



m 2 H ± 



-26, 

2^o 0*1 ~ A*3>> 

2«o()"2 - M3), 

r]oo + 



(2.29) 
(2.30) 
(2.31) 

(2.32) 



with 



m 2 h/ 



> m\„ 



(2.33) 



due to (2.16). The scalar sector also yields different interactions between the Higgs 



particles and gauge bosons. 



In the Yukawa sector of the MCPM (2.22) only the fermions of the third generation 



couple to cpi, and only tp\ has a non-zero vacuum expectation value. Hence only the 
t, b and r become massive through electroweak symmetry breaking, while the fermions 
of the first and second generation stay massless. The tree-level fermion masses in the 
MCPM are [I] 



m t 
m b 
m T 



vqc u 

vpci 
V2' 



mj = for all other fermions /. 



(2.34) 

(2.35) 

(2.36) 
(2.37) 



Similar to the SM Higgs boson, the p' couples to all massive fermions, i. e. to the t, b 
and r. However, the other Higgs particles behave differently. In (2.22) the second scalar 
doublet if2 couples only to the fermions of the second generation. Therefore its quantised 
excitations, the h\ h" and H^, couple to the c, s and //, although these particles are 
massless in the MCPM. The Feynman diagrams of the Higgs-fermion couplings can be 



found in fig. 2.1 
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2.3. The MCPM with three generations 




C,8,H 



ff ti, h" H~ 

Figure 2.1.: Couplings of Higgs bosons to fermions in the MCPM. While the p' couples to 
the third generation, the other Higgs bosons only interact with the massless 
second-generation fermions. 



The strength of these interactions is even more surprising. Because the coupling of 
ipi to the third generation fermions and the coupling of q>2 to the second generation 



fermions in (2.22) share the same coupling constants q, c u and c<i, the corresponding 
interaction vertices are weighted with the third-generation fermion masses. This is very 
different from the SM, where the second-generation fermions couple to the Higgs boson 
proportionally to the second-generation fermion masses, i. e. a lot more weakly. 

Finally, the CKM matrix turns out to be the unity matrix [T] , the quark mass eigenstates 
are weak interaction eigenstates as well. This again is not like in the SM, where the CKM 
matrix can be arbitrary (it is parametrised by four free parameters). 



2.3. The MCPM with three generations 



2.3.1. Particle properties and couplings 



Now the phenomenological properties of the MCPM can be put together. A list of all 
the particles and their key characteristics can be found in appendix A. 2 For the most 



important Feynman rules see appendix |A.3 There are four key differences between the 
MCPM and the SM: 



1. In the MCPM, the first two generations of fermions are massless at tree level. In 
the SM all fermions except the neutrinos are massive. 

2. In the MCPM, the CKM matrix is restricted to be unity. In the SM it is an 
arbitrary unitary matrix. 

3. There are five physical Higgs bosons (two of which are charged) in the MCPM 
compared to one in the SM. 
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2. The maximally CP-symmetric two-Higgs-doublet model (MCPM) 



4. Four of these Higgs bosons couple to second-generation fermions, but weighted 
with the third-generation fermion masses. There is no equivalent structure in the 
SM. 

Looking at the first two points, the MCPM does not seem to describe nature very well. 
However, the theory gives quite a good approximation. While there is no doubt that 
there are massive fermions in the first two generations, they still are a lot lighter than 



the t, b and r, see (1.1). The CKM matrix is certainly not quite 1, but close to it, see 
(1.2). So the MCPM approximates nature not too badly with relatively few parameters 
(11 compared to 18 in the SM). Generalised CP transformations might be approximate 
symmetries of nature. They might even be exact symmetries that are broken in some 
way. Finally there is also the possibility that radiative corrections generate fermion 
masses for the first two generations. Thus although there are some obvious spots of 
bother, this theory should not be discarded too easily and can provide a starting point 
for more research on the CP transformation behaviour of our world. 

So far the last two points of the list given above have neither been observed nor excluded 
by experiment. To discover or discard the MCPM as an exact, broken or approximate 
symmetry, they have to be studied carefully and checked at collider experiments like the 
LHC. 



2.3.2. Parameter space 



There are 11 coefficients in the MCPM Lagrangian before symmetry breaking, not count- 
ing the representation assignments: 

• Gauge couplings: g, g', g s 

• Higgs potential: m, ji 2 , (J>3, £o, Voo 

• Yukawa coefficients: c u , c<2, Q 



With (2.34) to (2.37) and (2.29) to (2.32), these and all other relevant variables such 



as vq can be expressed through measurable quantities. A useful set consists of particle 
masses, the Fermi constant Gp and the electromagnetic and strong coupling constants 
a and a s . The list of independent parameters of the MCPM then reads: 

• Couplings and gauge bosons: a, Gf, a s , mz 

• Fermion masses: mt, rribi m r 

• Higgs masses: m p r, m^i, mh", m H ± 
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2.4. Phenomenology at the LHC 



See appendix A. 4 for a list of important quantities expressed through these independent 
parameters. 

The gauge field properties and the fermion masses have already been measured exten- 
sively and with high precision, see [2] for a summary. While there is no direct evidence 
for any of the Higgs particles so far, electroweak precision measurements can be used to 
exclude a part of the parameter space for the Higgs particle masses [10 . In addition, 
results from the search for the SM Higgs boson probably apply to the p' as well, yielding 
strong restrictions on m p i. 



2.4. Phenomenology at the LHC 



2.4.1. Higgs production 



In proton-proton collisions, e. g. at the LHC, Higgs bosons can be produced via different 
processes. The most obvious one is quark-antiquark fusion, which is very similar to 
the Drell-Yan process [14, section 17.4]. The corresponding Feynman diagram is shown 



in fig. 2J2 As discussed in section |2.3| the p' couples only to t, b and r, while the 
other Higgs bosons couple only to the second-generation fermions. Hence the following 
quark-level processes are possible (if energetically allowed) : 

t + i-^p', c + c-^h', c + c^h", c + s^H + , 

b + b^-p, s + s-^/i', s + s^-h", s + c^-H^. 



P X 




P X 



Figure 2.2.: Production of Higgs bosons via quark-antiquark fusion. X denotes the rem- 
nants of a proton, experimentally seen as jets. 
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In [7, section 4.1] the cross-sections for these kind of processes have been calculated as 
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Here ^/s is the center-of-mass energy and 



w f 2 , 2 




dx 2 N q (xi)Nq>(x2)5 X\X2 



(2.38) 
(2.39) 
(2.40) 
(2.41) 

(2.42) 



where N q (x) is the parton distribution function (PDF) of the proton for quark q. 



Other production mechanisms include gluon-gluon fusion, see fig. 2.3 For the SM Higgs 
boson this is the dominant production process |16| . The p' behaves in the same way. 
However, it turns out that for all the other MCPM Higgs bosons quark-antiquark fusion 
is dominant |7J. This difference in behaviour is caused by the fact that the h' , h" and 
couple to the second-generation fermions, which are not that much suppressed in the 
PDFs, but the coupling strength is proportional to the heavy third- generation fermions. 
Hence quark-antiquark fusion is a lot stronger for these Higgs bosons than for the SM 
Higgs boson or the p'. 




Still, the production cross-sections are rather small. Evaluating (2.38) to (2.41) at the 
LHC design energy ^/s = 14 TeV and for Higgs masses around 150 GeV yields cross- 
sections for each Higgs boson with an order of magnitude of 4000 pb (where the 
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are slightly more likely to be produced than the ti and h"). For heavier Higgs bosons 
around 400 GeV this value drops to mere 80 pb. For exact numbers see fig. 12] or 



appendix C.l 



2.4.2. Higgs decays 



The Higgs bosons can decay in a number of ways. Similar to Higgs boson production, the 
most important modes stem from the fermion-fermion-Higgs vertices^] [7 section 3.1]: 



p - 


-> t + t, 


ti - 


■> c + c, 


/l"- 


■> c + c, 




p - 


+ b + b, 


ti - 


■> s + s, 


/l"- 


4 s + s, 




1 

p - 


+ r~ +r + , 


ti - 




/l" - 


4 M~ + M + > 





The corresponding Feynman diagram is given in fig. 2.4 

f 



c + s, 
s + c, 
M + + ^> 
M~ + *V 



/ 

Figure 2.4.: Higgs boson decay into fermion-antifermion pair. 
According to j7] , the decay width of such a process is given by 
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where the vertex coefficients a, b and the colour factor N c depend on the process and 



|6| S 



(2.43) 



are given in table A. 6 in appendix A. 5 These factors are proportional to the third- 
generation fermion masses. Because the t is significantly heavier than the b and r, the 
t and c channels dominate the decays, while the down-type quark and leptonic modes 



6 Of course these decays are only possible if the Higgs bosons are heavy enough. Given current restric- 
tions on the SM Higgs boson mass and taking into account that the p' behaves in a very similar way, 
this is probably not true for the p' — > t + i mode. 
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are suppressed by factors of roughly 10 3 to 10 5 . The total decay width amounts to 
approximately 9 GeV for a Higgs mass of 150 GeV or 24 GeV for heavy Higgs bosons of 
400 GeV. 

In addition to these simple fermionic channels, Higgs bosons can decay into other Higgs 
bosons and gauge bosons. In [TJ section 3.2 to 3.4] it has been shown that these decays 
are irrelevant for light Higgs bosons, but become important for Higgs boson masses 
above roughly 400 GeV. However, analysing and simulating these channels goes beyond 
the scope of this project. 



2.4.3. Angular distribution and resonances 

Apart from total cross-sections, the dependency on angular variables such as the polar 
angle #, the pseudorapidity rj or the azimuthal angle <p as well as the invariant mass of 
the decay products of the Higgs bosons is interesting. 

The Higgs bosons in the MCPM are scalars. In their rest frame, there is no preferred 
spatial direction and their decay is isotropic. For reactions involving production and 
decay of Higgs bosons this implies 

dcr da 

- — = const, — = const (2.44) 

dtp* d cos 0* 

where (p* and 9* denote the angles in the rest frame of the virtual Higgs particle. 

Note that this need not hold true in the lab frame. The quarks forming the Higgs 
bosons can each have arbitrary momenta^] so the virtual Higgs particle need not be at 
rest in the lab frame. In fact, it is usually moving fast in either of the beam directions. 
Therefore the angular distribution of its decay products is distorted and peaks in the 
beam directions. 

In processes including production and decay of an unstable particle such as a Higgs boson, 
the intermediate particle is not necessarily on-shell. Using the energies and momenta of 
its decay products, the invariant mass squared 

p 2 = E 2 - p 2 (2.45) 

of a virtual Higgs boson can be reconstructed easily. The dependency of the cross- 
section on this invariant mass is described by the relativistic Breit-Wigner distribution 
PH section 2.4.4] 

da k . , , 

(2.46) 



dy/p* (p 2 — m 2 ) 2 + m 2 T 2 
7 Up to an upper limit. This is described by the parton density functions. 
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Here m and T are the mass and decay width of the Higgs boson involved in the process, 
and A; is a constant of proportionality. This function peaks at p 2 = m 2 , it describes the 
resonance shape of the process. Note that this is only an approximation valid in the 
limit of narrow widths, r <C m. For broad peaks such as the Higgs processes of interest 
there will be some deviations. Most notably, the width of a broad resonance peak will 
slightly differ from the decay width of an on-shell intermediate particle [see for instance 
HU section 7.3]. 



2.4.4. Standard model background 

The MCPM Higgs bosons will never be produced on their own without any background; 
there are a lot of ordinary QCD and QED processes with the same final state that 
compete with Higgs production and decay. 

For the leptonic channels 



p + p 
p + p 



T~ +T++ X, 



the most important SM background is the Drell-Yan process consisting of annihilation 
of a quark-antiquark pair into an excited photon or Z boson, which decays into a pair 



of leptons |14| section 17.4]. The corresponding diagram is shown in fig. 2.5 There are 



similar diagrams for the /U~P^ and \i channels involving virtual W bosons. All of 
these processes have large cross-sections. However, leptons can be detected efficiently, 
and it is easy to reconstruct the four-momentum and thus the invariant mass of the 
virtual particle in such a process. 




Figure 2.5.: Drell-Yan process. Quark and antiquark from colliding protons annihilate 
into a virtual photon or Z boson, which decays into a pair of leptons. 



The quark channels are more complex. In addition to Drell-Yan-type processes, there 
are various QCD processes that yield two jets, see for instance fig. 2.6 Cross-sections 



of the strong interaction are typically large. The situation gets even worse when taking 
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hadronisation, jet tagging and reconstruction into account, which are complicated and 
therefore difficult to simulate and challenging experimentally. 



P X 




P X 
Figure 2.6.: Example for QCD background. 

So while the Higgs bosons of the MCPM decay mostly into quarks and much less into 
leptons, the huge background and the experimental difficulties of hadronic channels 
suggest that the leptonic channels are more promising [7]. Of these, the muonic events 
are most interesting, as the h', h" and couple to them in a way that is unique to the 
MCPM. 

While there is certainly a large background hiding possible MCPM signals, the situation 
is not hopeless. The gauge bosons involved in the Drell-Yan-type processes have spin 1, 
while the Higgs bosons are scalars. This leads to a different angular distribution of the 
leptons or jets. In their rest frame, virtual photons and Z bosons decaying to leptons 
approximately lead to a (1 + cos 2 t?*) distribution [21 section 25.4], which peaks in the 
beam directions. On the other hand, spin-0 particles decay isotropically as explained 
above. So by only looking at the central region of the detector, the SM background is 
reduced more strongly than the MCPM signal events. 




P 



Figure 2.7.: The full diagram describing production and muonic decay of the h! and h" . 
There are similar diagrams involving H^. 

In addition to separating the MCPM signature from the SM background, one also has 
to differentiate the MCPM signals from other SM extensions. Here the characteristic 
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generation structure of the MCPM is important: While the \x /j> + channels are relatively 
strong, the Higgs decays into other leptons are suppressed. 

Taking everything into account, the MCPM has a unique phenomenology. The p' behaves 
in a similar way as the SM Higgs boson and is therefore of less interest. For the h 1 , h" 
and H^, the main production channel at the LHC is predicted to be quark-antiquark 
fusion. These Higgs bosons decay mostly into c and c quarks, but the resulting jets are 
hard to tag and reconstruct and there is a huge QCD background. Therefore the most 
promising channels are the muonic Higgs decays as shown in fig. |2.7[ where the main 
background comes from Drell-Yan processes involving virtual photons or Z bosons. 
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generation 



It is possible to calculate the production cross-section for Higgs bosons at the LHC 



using (2.38) to (2.41). With (2.43) and similar expressions for other decay modes, the 
branching ratios for the different Higgs particles can be calculated. Assuming that the 
Higgs bosons are on-shell (i.e. "real" particles as opposed to virtual ones), the cross- 
section for the full process p + p — > H + X — > f + f + X is given by 

f full process = /^H ' -6// (3.1) 
H 

where an is the production cross-section for Higgs boson H and B^j is the branching 
ratio of the decay H — > f + /. This is the narrow- width approximation and usually 
valid if the intermediate particle is much heavier than the initial and final particles 
|18| . This assumption is likely justified for the MCPM, which allows the calculation of 
cross-sections for the MCPM without any Monte-Carlo simulation, as has been done in 

m- 

However in such a calculation a lot of effects cannot be easily taken into account, such as 
hadronisation and detector behaviour. For a full analysis a different approach is needed: 
the random generation of events with Monte-Carlo methods followed by the simulation 
of QCD and detector effects. This approach yields data that is directly comparable to 
experimental results. 

The main part of this project was the implementation of such a Monte-Carlo event 
generation for the MCPM. The framework of MadGraph 5 has been used, which will 
be explained below. Afterwards the key features of the MCPM implementation are 
described. 



3.1. Introduction to MadGraph 5 

MadGraph 5 |11] is the latest version of the matrix element generator MadGraph. It is 
bundled with the event generation package MadEvent. For simplicity, both will be re- 
ferred to as MadGraph in the following. Given an arbitrary process for any implemented 
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model, the program lists all contributing tree-level diagrams. For each of these diagrams, 
events are randomly produced and the cross-section or decay width is calculated. 

The output of such a simulation is a list of events in the Les Houches Event File (LHEF) 
standard [E]. For each event, it includes the particles involved and their four-momenta. 
This can be read by programs such as PYTHIA [20 , with which hadronisation effects 
can be simulated. It is also possible to perform a detector simulation on the results, 
for instance using GEANT [21]. Event generation with MadGraph can thus be used as 
the first part of a complete Monte-Carlo analysis leading to a comparison of theory and 
experiment — and ultimately the discovery or exclusion of new physics models such as 
the MCPM. 

MadGraph is a very flexible framework and allows the implementation of new models in a 
simple way using the Universal FeynRules Output (UFO) format [22] . Its major downside 
is the limitation to tree-level processes. However as seen in chapter |2.4| the interesting 
and relevant processes do not involve any loops, so this drawback is acceptable. 



3.2. MCPM implementation 

In the course of this project, I have implemented the MCPM in MadGraph using the 
UFO interface [22]- Basically, the model is described by three lists: 

• A list of the parameters of the MCPM. As independent parameters, the quantities 



listed in section 2.3.2 were used. For these variables, default values were cho- 
sen according to PDG recommendations [2]. For the Higgs boson masses, values 
consistent with the restrictions presented in [10J were set. 

• A list of the particles of the MCPM and their properties, including spin, mass, 
electromagnetic charge and colour charge. 

• A list of the interaction vertices in the MCPM. For each vertex, the participating 
particles and the corresponding Feynman rule had to be given. Some important 
Feynman rules have already been calculated in [7], others had to be derived first. 

The MCPM model files are written in Python and span roughly 1300 lines. For more 



details, see appendix B.2 



All particles and couplings of the MCPM are included. This not only allows the analysis 
of signal processes involving Higgs bosons, but also the simulation of any tree-level 
background process. There is one important feature: strictly speaking, the MCPM 
forbids masses for the first two generations of fermions at tree-level. However nature 
seems to disagree with this prediction. This suggests that generalised CP transformations 



might only play the role of approximate symmetries, as discussed in section 2.3. It may 
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be sensible to analyse a model that combines the scalar sector of the MCPM with the 
experimental evidence for fermion masses and the CKM matrix [7]. For this reason the 
MadGraph implementation of the MCPM lets the user set fermion masses for the first 
two families as well. This affects phase space calculations, but the interactions remain 
unchanged. In a similar way the CKM matrix can be set to an arbitrary matrix, which 
changes SM processes, but not the Higgs boson vertices. With these parameters the user 
can decide whether to analyse the strict MCPM as described in chapter [2] or a hybrid 
theory with fermion masses and a SM-like CKM matrix. If any SM processes play a role 
(for instance for background simulation), the latter seems to make more sense and will 
be used in the following. 

Now the model can be used in the same way as other models are used in MadGraph. 
In order to generate events for a process, only a few commands are necessary. As an 
example, the lines 



import model MCPM -modelname 
generate p p > hi > m- m+ 
output -f 
launch 



lead to the generation of events for the process p + p — > h' + X — > \x + p + + X. The 
free parameters can be set in a configuration file. For detailed instructions on how to 



use the model, see appendix B.l 



30 



4. Validation 



The MadGraph implementation of the MCPM was validated in two different ways. First, 
total cross-sections and decay widths were checked for a number of different scenarios. 
For this purpose Higgs-boson production and decay were treated independently. Second, 
angular dependencies and invariant mass distributions were analysed for full processes. 



4.1. Production cross-sections and decay widths 
4.1.1. Methodology 

The Monte-Carlo simulation was validated by simulating different processes with Mad- 
Graph and comparing the resulting cross-sections as well as decay widths to the corre- 
sponding theoretical values. Production and decay were treated independently, i. e. the 
Higgs bosons were assumed to be on-shell. 



As discussed in section 2.4 the relevant production processes are 

p + p -5- ti + X, 
p + p -> h" + X, 
p + p^ H~ +X 

via quark-antiquark fusion. As before, X denotes the fragmented remnants of the pro- 
tons. The interesting decay modes include 

h — > c + c, h — > c + c, H~ — > s + c, 

These three production and six decay processes were simulated. 

For the MCPM parameters, PDG recommendations [2] were used. Of course this leaves 
open the Higgs masses. Here four different points in parameter space have been chosen 
that fulfill two conditions: 

• They cover a wide range of the energies relevant at the LHC. In addition the ratios 
of the Higgs masses to each other are varied. 
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All parameter sets agree with the electroweak precision measurements discussed 
in [10] within one standard deviation. They also agree with the stability criterium 
d2l33b. 



The chosen Higgs mass values are given in table 4.1 



Parameter set m p > [GeV] [GeV] mh" [GeV] m H ± [GeV] 



1 


125 


100 


100 


150 


2 


125 


400 


200 


250 


3 


170 


200 


150 


200 


4 


170 


300 


250 


300 



Table 4.1.: Parameter sets used for the validation. 



The center-of-mass energy was set to y/s = 14 TeV and no phase-space cuts have been 
applied. For the decays, 20,000 events have been produced in each MadGraph run. The 
production processes are more complex as they include multiple diagrams weighted by 
the PDFs of the proton. Here MadGraph was set to simulate 100,000 events in each run. 
The PDF CTEQ5L (23] was chosen for compatibility reasons. 

The theoretical values were calculated for the same processes and settings. While the 



calculation of the decay widths (2.43) was straightforward and was done analytically, 



the production cross-sections (2.38) to (2.41) include the PDF and hence had to be 



calculated numerically, using the CTEQ5L Fortran libraries 
For each of the nine processes given above and each of the four parameter sets in table 



4.1 , a cross-section or decay width was simulated with MadGraph and a theoretical value 
was calculated. So in total there are 36 data points, each allowing the comparison of the 
Monte-Carlo result with the expectation. 



4.1.2. Results 



The results of this procedure can be found in appendix C.l Fig. 4.1 shows residual 



plots, i.e. the deviations of MadGraph and theoretical results. Without a doubt, the 
Monte-Carlo results agree well with the expectation. For the decays, the numbers are 
exactly the same in all the five digits MadGraph indicates. The largest deviation in the 
production cross-section data is only 0.06 per cent, which is still well below one standard 
deviation of the Monte-Carlo results^] The MCPM model in MadGraph passes its first 



1 This is even better than one might naively expect. Should there not be some statistical fluctuations 
leaving only 68% of data points within one standard deviation of the theoretical values? The answer 
lies in the internal mechanisms of MadGraph The software calculates the cross-sections for 
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test. 



4.2. Event distributions 

After validating the MadGraph results on the level of production cross-sections and 
decay widths, full processes including virtual Higgs bosons were simulated. The angular 
distributions and resonance shapes of the decay products were analysed. 



4.2.1. Methodology 

The processes 

p + p -+ c + c + X, (4.1) 

p + p^ s + c + X, (4.2) 

p + p -+ pT + /i + + X, (4.3) 

p + p-> fT + Up + X (4.4) 



of the form shown in fig. 2.7 were simulated with MadGraph. Only MCPM Higgs bosons 



were allowed as intermediate particles, no SM background was included. 



The parameters were chosen as in section |4~T| with three exceptions: every run included 

100,000 events, and the newer PDF CTEQ6L [21] was used. The Higgs masses were set 
to 

mf, = 125 GeV, m# = 200 GeV, (4.5) 

m£„ = 150 GeV, m^± = 150 GeV. (4.6) 

This parameter set will be referred to as set A in the following. It complies with the 
restrictions found in |10| . 

For these processes and parameters, MadGraph calculates the following total cross- 
sections: 

a c5 = 4980 pb, (4.7) 

a s - c = 4760 pb, (4.8) 

a^-^ = 0.175 pb, (4.9) 

= 0.168 pb (4.10) 



individual quark-level processes precisely, the only uncertainty is introduced through the PDFs. 
However it gives out a standard deviation that does not represent this uncertainty of the calculated 
cross-section or decay width, but the magnitude of statistical fluctuations that are to be expected 
for an experiment observing the number of events set in MadGraph. 
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(a) Higgs boson production 
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(b) Higgs boson decay 



Figure 4.1.: Residual plots showing the deviations of the Mad Graph results and theoret- 
ical expectations. Good agreement is found for each process and parameter 
setting. 
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with negligible statistical errors. 

In addition to calculating cross-sections, the resulting event files were examined using the 
tool MadAnalysis, which is included in the MadGraph 5 package. It allows the analysis 
of arbitrary distributions based on the four-momenta of the particles. In the following, 
angular dependencies and invariant mass distributions are shown. As MadAnalysis did 
not provide all the features needed, some routines had to be added by hand, including 
the correct handling of statistical uncertainties. ROOT [25] and GNUPlot |26j were used 
for fitting and plotting. 



4.2.2. Angular distribution 



As described in section 2.4.3 the Higgs bosons should decay isotropically in their rest 
frame, which is the same as the center-of-mass frame of the decay products. The dis- 
tribution is different in the lab frame. For the four Monte-Carlo samples generated, the 
angular distributions in both reference frames were analysed, some of the results are 
shown in fig. 4.2 and fig. 4.3. More plots can be found in appendix |C. 2 



For the cost?* and (p* distributions in the Higgs boson rest frame, constant-value fits 
were performed. They agree well with the data, the \ 2 values are listed in table 4.2 



Channel 


cos-d* distribution 
X 2 1 ndf 


tp* distribution 

X 2 1 ndf 


cc 


40.7/49 


53.8/61 


sc. 


53.3/49 


58.2/61 




59.9/49 


79.1/61 




53.3/49 


58.2/61 



Table 4.2.: Constant fits to the angular distributions in the Higgs boson rest frame. 



4.2.3. Resonance shapes 



In addition to angular distributions, the dependency of the cross-section on the invariant 
mass of the Higgs boson is important. For each possible intermediate particle a reso- 
nance peak at its mass is expected. Its form should approximately be described by the 



relativistic Breit-Wigner distribution (2.46), see section 2.4.3 



For the four channels given above the invariant mass distribution was analysed. Fig. 4.4 
shows the resonance shape in the channel, the others are similar. Breit-Wigner 

distributions were fitted to the data with good agreement. The resulting fit parameters 
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were then compared to the Higgs boson masses (4.6) and the Higgs boson decay widths 



(2.43 ), which were calculated for the dominant fermionic decay modes (see section 2.4.2 ) 



The results can be found in table 14.3 



Channel 


m [GeV" 
Fit 


Theory 


T [GeV] 
Fit 


Theory 


X 2 1 ndf 


cc 


149.90 ±0.02 
199.87 ±0.04 


150.00 
200.00 


9.00 ±0.04 
11.24 ±0.08 


8.83 
11.77 


515.5/294 


sc 


149.94 ±0.02 


150.00 


8.81 ±0.03 


8.83 


726.3 / 297 


U.-U.+ 


149.90 ±0.02 
199.85 ±0.04 


150.00 
200.00 


9.02 ±0.04 
11.29 ±0.08 


8.83 
11.77 


464.6 / 294 




149.91 ±0.01 


150.00 


8.76 ±0.03 


8.83 


699.5 / 297 



Table 4.3.: Breit-Wigner fits to the resonance shapes. 

There are some small deviations. The reduced \ 2 values are in the region of 2, and 
the fitted parameters slightly differ from the Higgs boson masses and decay widths. As 
explained in section 2.4.3| this is not surprising and a consequence of the relatively large 



widths of the Higgs bosong^j 

Taking everything into account, the MCPM implementation produces exactly the ex- 
pected angular distributions and resonance peaks. Together with the cross-section check 



in section 4.1 this result concludes its validation. Due to time constraints, not every 
process and not every parameter configuration have been tested. Nevertheless, I consider 
the model to be validated satisfactorily in order to be used for analysis. 



2 Considering that the bin size used for this analysis is 0.5 GeV, binning effects also seem possible. 
However decreasing the bin size does not change the results by a relevant amount. 
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(a) Distribution in the lab frame. 
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(b) Distribution in the center-of-mass frame of the muons, which is the rest frame 
of the intermediate Higgs boson. 

Figure 4.2.: Dependency of the cross-section on the polar direction •& of the \x~ . 

•d is defined to satisfy cos 1 !? = ±1 in beam direction, cosi? = represents 
a muon moving perpendicular to the beam axis. In the lab frame, the 
distribution is strongly peaked in beam direction due to the momentum of 
the virtual Higgs boson. After a boost into the center-of-mass frame of the 
muon pair, the distribution becomes isotropic. 
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Figure 4.3.: Dependency of the fj,~ /U + cross-section on the azimuthal direction cp of the 
pT . To leading order, this distribution is the same in the lab frame and in 
the center-of-mass frame of the muons and should be constant. The Monte- 
Carlo data agrees very well with this expectation. 
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Figure 4.4.: Dependency of the n. fi + cross-section on the invariant mass of the muon 



pair. This is identical to the invariant mass (2.45) of the virtual Higgs 



boson. There are resonances of the h! at 200 GeV and of the h" at 150 GeV. 



A superposition of two Breit-Wigner shapes (2.46) is fitted to the data with 
good agreement. 
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5. Analysis of MCPM signatures 



The MCPM MadGraph implementation can now be used to analyse an important ques- 
tion: can the MCPM be discovered or excluded at the LHC? A full answer requires 
a thorough analysis including the simulation of fragmentation and detector behaviour, 
which goes well beyond the scope of this thesis. For the hadronic channels any su- 
perficial analysis using only quark-level events is meaningless, as hadronisation and jet 
reconstruction are complicated and important. These channels will not be further anal- 
ysed. 

However there is no hadronisation for leptons, and muons can be reconstructed relatively 
precisely. For the /i~ /U + channel the hard-process events generated by MadGraph are 
worth a look even without fragmentation and detector simulation. Using the new Mad- 
Graph implementation, the MCPM signatures are first compared to the SM background, 
both without and with phase-space cuts. Then ATLAS results from 2011 are checked 
against the MCPM predictions. 



5.1. Comparison of MCPM signals to SM background 

In this section the MCPM signal events are compared to the SM background for the 
/U~ fi + channel. First the raw signals are analysed. Then basic r/ and px cuts motivated 
by detector properties are applied. 



5.1.1. Raw signal 



For this analysis, both signal and background events were generated with MadGrapfj^] 
The relevant signal processes are once again the production of h' and h" bosons via quark- 
antiquark fusion with subsequent decay into fj> + pairs, see fig. 2.7 As described in 



section 2.4.4 the background is dominated by the 7* and Z Drell-Yan processes. These 
lead to the same final state as the signal procesess and thus cannot be reduced easily. 



1 Signal and background can be treated separately. Because of the different chirality structures, there 
is no interference between signal and background diagrams. 
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5. Analysis of MCPM signatures 



The same parameters as in section 4.2 were used. This includes the Higgs boson mass 
set A (4.6) with mfy = 200 GeV and m^„ = 150 GeV. In order to analyse a broader 
spectrum of resonances, a second Higgs mass set B was chosen in addition: 



m B p , = 125 GeV, 



m 



B 



m 



B 

h" 



250 GeV, 



in 



B 



400 GeV, 
300 GeV. 



(5.1) 
(5.2) 



For the background simulation, more than 500,000 events were generated in each run. 
All simulations were performed for the LHC design energy y/s = 14 TeV as well as its 
current center-of-mass energy y/s = 7 TeV. In the following the plots for 14 TeV are 
shown, the corresponding 7 TeV plots can be found in appendix |D.l 



Fig. 5.1a shows the signal and background cross-sections. The resonance peaks of the 
Higgs bosons are tiny in comparison to the background. This can be quantified by 
comparing the MCPM and SM cross-sections for the invariant mass rangd^] between 
run — r# and run + One finds 



""signal 
""background 



resonance 




for yfs 
for yfs 



14 TeV 
7 TeV 



(5.3) 



if the Higgs bosons are in the range of 150 to 250 GeV. For heavier Higgs bosons these 
factors decrease. See appendix |D.l for more detail. 



5.1.2. Imposing cuts 

There is no detector that is able to detect and reconstruct every muon. For instance, in 
the ATLAS detector at CERN the muon trigger chambers only cover the pseudorapidity 
region \rj\ < 2.4 [27\ chapter 6]. In addition, muons are efficiently reconstructed only for 
high transverse momenta pr- It is therefore sensible to impose cuts on the Monte-Carlo 
events. Starting from the event samples used in the last section and requiring 

\n\ < 2.4, (5.4) 
p T > 15 GeV (5.5) 

should give a more realistic picture of MCPM discovery potential. 

While this of course decreases statistics, there is one advantage. As discussed in section 



2.4.4 the spin-0 Higgs bosons are more likely to decay into this central region than the 



2 By choosing a smaller invariant mass bin, e.g. mu ± \Vu, the ratio — tT °' 8 " a ' — increases. But from 

^ ^background 

an experimental point of view it makes sense to use relatively large mass bins in order to get more 
statistics and due to detector resolution. 
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5.1. Comparison of MCPM signals to SM background 



spin-1 Z bosons and photons which make up the background. Hence these cuts reduce 
the background more drastically than the signal events. 



Fig. 5.1b shows the resulting invariant mass distribution. As expected, the MCPM 
signature is more clearly visible. For Higgs masses around 150 to 200 GeV and looking 
only at the resonance regions where the invariant mass of the muons is in the range of 
m H ~ r# to m H + T H , 

"signal 



""background 




s = 14 TeV , , 

(5.6) 

resonance I 3 • 10 A for J~S = 7 TeV. 



Again, this value decreases for heavier Higgs bosons. More data can be found in appendix 

ED 

Now one can calculate how much statistics is needed so that the MCPM resonance peaks 
are locally significant compared to the background, i. e. larger than na, where a is the 
statistical uncertainty of the background. Conventionally n = 5 is considered a discovery, 
while n = 2 is often used for exclusion limits^} I calculated the integrated luminosities 



/ Ldt needed for local 2a and 5a significance. The results are shown in table 5.1 



m H [GeV] 




/ Ldt [fV 


- 1 ] 






2a significance 




5a significance 






y/s = 7 TeV ^s = 


14 TeV 


y/s = 7 TeV y/s = 


14 TeV 


150 


6.7 


1.2 


41.9 


7.7 


200 


22.0 


3.1 


137.7 


19.1 


250 


68.5 


7.9 


427.9 


49.4 


400 


866.4 


85.9 


5415.0 


536.8 



Table 5.1.: Integrated luminosities that are needed for local significances of the Higgs 
boson resonances. 

I have to stress though that these luminosities are very rough estimates. Only local 
significance in one channel is described, the combination of different peaks and different 
channels will improve the situation a bit. On the other hand, these numbers do not in- 
clude detector efficiencies and resolution, only background and signal tree-level processes 
are taken into account. The experimental reality is much more complicated. 

Keeping this caveat in mind, these numbers may nevertheless hint at the order of mag- 
nitude of statistics needed for MCPM signatures to become visible. As of October 31, 

3 These values are valid only locally, i. e. in the invariant mass range corresponding to the resonance 
peaks, and they cannot be simply translated into exclusion limits based on confidence levels. More 
sophisticated statistical tools are needed for this. 
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5. Analysis of MCPM signatures 



2011, the ATLAS detector at CERN collected a total of 5.25 fbr 1 of data at y/s = 7 
TeV [28]. It is expected that the LHC will collect significantly more data in 2012 and 
will reach its design energy of y/s = 14 TeV a few years later. Looking at the critical lu- 
minosities in table [STTj the next years will be an exciting time for the MCPM. Although 
a formal 5u discovery seems very unlikely to happen soon, there might be either hints 
at the MCPM or exclusion limits on its Higgs boson mass space. 



5.2. A brief look at the Z' search 

The MCPM is not the only extension of the SM that predicts fi~ u. + events, and due to 
reconstruction efficiency and low background this channel is used in various searches for 
signs of new physics, for instance Z' boson resonances (for an introduction to Z' physics, 
see [22]). In [3D], 1.21 fb -1 of data taken with the ATLAS detector at y/s = 7 TeV were 
analysed for excesses over the SM predictions. In contrast to the estimates in the last 
section, the Z' search includes a proper analysis of hadronisation, QCD background 
faking muons and detector properties. 

In |30| , the expected and measured number of events passing certain cuts were compared 
for different invariant mass bins. The most relevant cut is on the transverse momentum 
Pt- Each muon has to satisfy 

p T > 25 GeV. (5.7) 

Due to detector geometry, 

\r]\ < 2.4 (5.8) 

is required. Other cuts limit the position of the primary vertex and the impact param- 
eter. 

Using MadGraph, the cross-sections for the usual signatures involving h' and h" 

bosons were simulated. The parameters and cuts were chosen such as to closely resemble 
the parameters of the Z' searclj^J For the Higgs boson masses, two different sets of 



parameters C and D were chosen, see table 5.2 Notably, parameter set C represents a 



best-case scenario with a large Higgs boson resonance peak at 120 GeV. The technical 



parameters were chosen as in section 5.1 



Table 5.3 shows the results of the Z' search and the MCPM predictions for the two 
settings. Even for the best-case scenario C the MCPM predicts less of an excess than 
the statistical and systematical uncertainties of data and background analysis. 



4 Some cuts could not be replicated in MadGraph, for instance restrictions on the primary vertex 
position. 
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5.2. A brief look at the Z' search 





Parameter set 


m pl [GeV] 


m h , [GeV] 


•nihil [GeV] 


m H ± [GeV] 


C 
D 


125 
125 


120 
250 


120 
150 


125 
150 



Table 5.2.: Parameter sets used for the comparison of Z' search results and MCPM 
predictions. 



fjj ji + invariant mass bin [GeV] 


110 


-200 


200 


-400 


ATLAS data according to [3D] 


5406 


± 74 


557 


±24 


SM prediction according to [30] 


5434 


±150 


571 


±23 


MCPM signal (parameter set C) 


42 


± 6 





± 


MCPM signal (parameter set D) 


19 


± 4 


3 


± 2 



Table 5.3.: Number of events measured at ATLAS compared to SM predictions and 
MCPM signatures. Here the invariant mass bins are relatively large, but 
even with smaller bins the MCPM signature is insignificant compared to the 
background uncertainties. 

A thorough analysis including detector simulation will further decrease the MCPM pre- 
diction. Thus it can be concluded that the data presented in [3D] is probably not sufficient 
to discover or exclude the MCPM for Higgs boson masses above 100 GeV. Of course, the 
combination with other channels and the analysis of the full data taken in 2011 might 
yield a different picture. 
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5. Analysis of MCPM signatures 



1000 
100 
10 
1 

0.1 
0.01 



Signal A, enlarged by factor 10 ■ 
Signal B, enlarged by factor 10 

Background ■ 



50 100 150 



(a) Raw data - no cuts applied. 



200 250 300 
m(fi-,fi+) [GeV] 



350 



400 



450 



ft 



100 
10 

1 

0.1 
0.01 



50 



Signal A, enlarged by factor 10 l I 
Signal B, enlarged by factor 10 

Background ^^^^ 



\ti\ < 2.4, p T > 15 GeV 



100 



150 



200 250 300 

m{fi-,fi + ) [GeV] 



350 



400 



450 



(b) Imposing cuts on rj and px- 



Figure 5.1.: MCPM events compared to the SM background in the channel as a 

function of the invariant mass of the muon pair, y/s = 14 TeV. While the 
Z resonance at 91.2 GeV is clearly visible, the Higgs resonances would be 
hard to notice without the magnification factor of 10. 
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6. Summary 



In this Bachelor thesis I implemented a two-Higgs-doublet model with maximal CP 
symmetry, the MCPM, into the Monte-Carlo event generation package MadGraph. 

The MCPM is an extension of the Standard Model that is based on the requirement of 
invariance under generalised CP transformations. It features a unique fermion structure, 
in which only the third-generation fermions are massive. Of course this does not describe 
nature precisely, but the model gives a good approximation of what has been observed. 
The theory also predicts five Higgs bosons, of which four couple only to the second- 
generation fermions. At colliders such as the LHC, these Higgs bosons will be produced 
mainly via quark-antiquark fusion. Their fj,~ (i + decay modes seem most promising for 
a discovery. For a proper search for MCPM signatures, a Monte-Carlo simulation is 
needed. 

In the course of this thesis, the MCPM was implemented into MadGraph, allowing the 
calculation of cross-sections and the random generation of events for arbitrary tree-level 
processes. The output is in the LHEF format and can be used for further analysis, for 
instance using PYTHIA and GEANT. Therefore this MadGraph implementation can 
be used as the starting point for a full and thorough Monte-Carlo simulation of the 
MCPM. 

The implementation was validated successfully with different methods. It was then used 
to compare the MCPM signatures to the SM background at LHC energies. This analysis 
was only done for the / u^ / Lt + channel and restricted to hard processes, so the results are 
only a rough approximation. It was found that even after applying appropriate cuts, the 
Higgs boson resonance peaks are small compared to the background. For a Higgs boson 
of 150 GeV (200 GeV), 6.7 ftr 1 (22 fbr 1 ) of data at a center-of-mass energy ^ = 7 
TeV have to be collected so that a local 2a significance becomes possible, not taking 
detector inefficiencies into account. It seems that the exclusion of a part of the Higgs 
mass parameter range or the discovery of small local excesses hinting at the MCPM 
might be within reach in the next years. Finally the MCPM predictions were compared 
to 2011 results of the Z' search at the ATLAS detector, again hinting that more data 
must be collected in order to discover or exclude the MCPM for Higgs masses of more 
than about 100 GeV. 
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Appendix A. 

Further details of the MCPM 



A.l. The THDM Lagrangian 



The Lagrangian of a general THDM is given in (2.8) with the separate sectors given 



in (2.3), (2.5) and (2.7). The MCPM is more specific, its Higgs potential and Yukawa 



sectors can be found in (2.15) and (2.22). In these terms, the following quantities 
appear: 

• u m , d m , e m and v m are fermionic fields representing up- type quarks, down- type 
quarks, charged leptons and lepton neutrinos, m is a generation index, three 
generations have been observed in nature. 



With the projection operators 
_ IT 75 

L/ R - 7, 



(A.l) 



the left-handed and right-handed chiral projections of these fields can be defined 
as *1> L /R = Pl/r4>- 

The left-handed particles are organised in SU(2) doublets, 



(V- 



and If 



(A.2) 



There are two complex scalar doublets 



(A.3) 



Gt, WJj and are the gauge fields of the SU(3)c X SU (2)iX(7(l)y gauge group. 
g s , g and g' are the corresponding coupling constants. 
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Appendix A. Further details of the MCPM 



• Their field strength tensors are given by 

G%v = d»Gi - d v G\ - g s f ijk G^G k v (A.4) 

= W - d u W; - g^Wffi (A.5) 

B^ u = d^B u - d v B^. (A.6) 

• The covariant derivatives are 

D^qT a = (fy + i9.I>GW + igTW^ + i^B^A qf (A.7) 

D»l? = (dp + igrw; - i 9 -B^ If (A.8) 

DpuT = (a„ + ig.L'G** + i^-B^ (A.9) 

D^ = {d^-ig'B,)e m (A.ll 



D»dT =[d, + ig s VG^ - i»-B^\ d% p (A.10) 



Dtfi = [ + igVWl + ) (A.12) 



• L l and T l are SU{2>) and SU{2) generators in the fundamental representation. 

• f l i k and e^ k are the structure constants of these representations. 

• There are Higgs potential and Yukawa coefficients. In the general THDM these 
are called ay, b Mmn , C l ™ n , C d mn and C\ mn . In the MCPM with three generations 
the coefficients are fix, fj,2, M3, £o> VoOj c u , Cd and q. 

• Colour and spinor indices are suppressed (with the exception of the quark colour 
indices a in the covariant derivative). 

• Feynman slash notation is used: 

U) ■= 7 m £>m- (A. 13) 

A. 2. Particle spectrum 

All the particles of the MCPM and their most important properties and couplings can 



be found in table A.l table A. 2 and table A. 3 
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A. 3. Feynman rules 



Gen. 



Fermion 



Mass 



Charge 



Couplings 



Gauge 



Higgs 



u 





2 

I 

3 


G, W ± , Z, 7 


d 





G, W ± , Z, 7 


e 





-1 


W ± , Z, 7 










W ± , Z 





c 





2 

! 

3 


G, W ± , 




h>, 


h", 


H± 




s 





G, W ± , 




h>, 


h", 


i?± 


2 


A* 





-1 


w ± , Z, 


7 


h', 


h", 


H± 












w ± , z 




h', 


h", 









m t 


2 

1 
3 


G, W ± , 


z, 


7 


P' 




b 


m b 


G, W ± , 


z, 


7 


p' 


3 


T 


m T 


-1 


w ± , z, 


7 




p' 




U T 








w ± , z 






p' 



Table A.I.: Fermions in the MCPM, their key properties and couplings. Unlike in the 
SM, the first two generation of fermions are massless. The third-generation 



fermion masses are given in (2.34) to (2.37) 

A. 3. Feynman rules 



In appendix A] the MCPM Lagrangian after symmetry breaking is presented and 
the Feynman rules for several vertices are derived. The most important vertices and the 



corresponding expressions are given in table A. 4 and table A. 5 below 



A. 4. Parameters and their relations 



As described in section 2.3.2 , the independent parameters used in the analysis are a, Gf, 
a s , mz, rn t , nib, TOr, m p >, nth', mh" and m H ±. All other quantities can be calculated 
from these parameters as follows J7J appendix A] : 



Co 
Voo 

Mi 



1 2 



1 



2v, 



2 {m 2 H ±+m 2 pl 



1 ( 2 



m 



H± i ' 



(A.14) 
(A.15) 

(A.16) 
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Appendix A. Further details of the MCPM 



Boson Mass Charge Couplings 









Fermions 


Gauj 






Hi! 


ygs 




G 








Quarks 


G 












W ± 


mw 


±1 


All fermions 


W ± , 


z, 


7 


ft, 


h>, 


h", i?± 


z 


m z 





All fermions 


W ± , 


z, 


7 


p', 


h>, 


h", H± 


1 








Quarks, charged leptons 


W ± , 


z, 


7 


h>, 


h" 





Table A. 2.: Gauge bosons in the MCPM. The gauge bosons and their couplings to the 
fermions are exactly the same as in the SM. 



Boson 



Mass 



Charge 



Couplings 



Fermions 



Gauge 



Higgs 



p' 


m p > 





t, 


6, 


T 


W ± , 


z 




all 


b! 


my 





c, 


s, 


/' 


W ± , 


z, 


7 


all 


h" 


mw 





c, 


s, 


P 


w ± , 


z, 


7 


all 


H± 


m H ± 


±1 


c, 


s, 


P, "ft 


w ± , 


z, 


7 


all 



Table A. 3.: Higgs bosons in the MCPM. While the p' is similar to the SM Higgs boson, 
the other Higgs particles have very distinctive properties. The Higgs boson 



masses are given in (2.29) to (2.32). Note that some of the couplings are 
four-point vertices. For instance, the uncharged b! and h" couple to the 7 
only via vertices involving and bosons as well. 



P2 
P3 

e 

sin 2 Ow 



m 2 w 




v = 2-iG F \ 



(A.17) 

(A.18) 
(A.19) 

(A.20) 

(A.21) 
(A.22) 
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A.4. Parameters and their relations 




Table A.4.: Feynman rules for the most relevant vertices in the MCPM, part 1. 
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Appendix A. Further details of the MCPM 



Vertex Corresponding expression 




Table A. 5.: Feynman rules for the most relevant vertices in the MCPM, part 2. There 
are many more, including couplings of gauge bosons to Higgs bosons and 
Higgs boson self-couplings. The arrows on the lines denote the flow of 
negative charge. 
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A.5. Higgs boson decay modes 



A. 5. Higgs boson decay modes 



In (2.43) the width of Higgs decays into fermion-antifermion pairs was given, including 



coefficients that are specific to the different decay modes. These coefficients have been 



calculated by Maniatis and Nachtmann table 2] and are listed in table A. 6 below. 



Higgs boson 


Decay mode 


a 


b 


N c 


P 


it 


m t 





3 




bb 


m b 





3 




T - T + 


m T 





1 


h! 


ce- 


-m t 





3 




ss 


-m b 





3 




fl-fi + 


—m T 





1 


h" 


ce- 





im t 


3 




ss 





—imb 


3 









—im T 


1 


H+ 


cs 


m t -m b 


m t +m b 


3 
1 


V2 
m T 

V2 


V2 
m T 

V2 


H~ 


sc 


m t -m b 

V2 
m T 

V2 


mt+mb 

V2 
m T 

V2 


3 
1 


Table A.6.: Hig 


gs boson fermionic decay 


modes and their coefficients. The decay 


widths 



for these modes are given in (2.43) 
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Appendix B. 



Technical aspects of the MadGraph 
implementation 

B.l. Instructions on usage 

The MCPM model in MadGraph is used in the same way as every MadGraph model, 
see [TT] for an overview. First it has to be loaded. After starting the MadGraph binary 
in a shell this can be done with the command 

import model MCPM -modelname 

where the option -modelname is needed to ensure the correct particle names. Then 
processes are defined and the event generation started: 

generate p p > hi > m- m+ 
add process p p > h2 > m- m+ 
output -f 
launch 

The generate and (optionally) add process commands are used to specify the processes 
MadGraph has to evaluate, while output and launch launch the event generation. The 
particle names used in the implementation are given in the following section. Before 
events are produced, MadGraph asks the user whether the standard parameters are to 
be used or if he wants to modify them. 

During the simulation, an in-browser status page keeps the user informed. In the end, 
MadGraph gives out the contributing Feynman diagrams, the cross- sections or decay 
widths, the corresponding uncertainty and a compressed file containing the generated 
events in LHEF [19 format. This can be used for further analysis, for instance with 
MadAnalysis or PYTHIA [2D]. 
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B.2. List of particles and parameters 



B.2. List of particles and parameters 

A list of MCPM particles and parameters with their MadGraph names are presented in 
table B.l and B.2| respectively. 



Particles 


Names in MadGraph MCPM implementation 


Protons 


P 


Quarks 
Antiquarks 
Leptons 
Antileptons 


u d c s t b 
u- d- c~ s~ t~ b~ 
e- m- tt- ve vm vt 
e+ m+ tt+ ve- vm- vt- 


Gauge bosons 


a z w+ w- g 


(/, h', h" 


rho hi h2 
h+ h- 


Table B.l.: The MCPM particles and their names in the MadGraph implementation. 



The parameters may be a bit confusing, because the c, e, [i can be set to be massive, 
while the u, d and s are always massless. This does not have anything to do with 



the MCPM. As explained in section 3.2 some of the consequences of strict symmetry 
under generalised CP transformations have been dropped in this implementation, so it 
is generally possible to set masses for the first- and second-generation fermions. How- 
ever, MadGraph treats the very light fermions as massless (even in the SM, where the 
theory certainly predicts something else). This does not change the physics at collider 
experiments in the TeV range, but it saves some computation time. Therefore the list 
of parameters includes the masses of some fermions of the first two generations. 

MadGraph also allows to set decay widths for all particles, which are named similarly 
to the mass parameters. This is important for further analysis with programs such as 
PYTHIA, but these parameters do not influence the generation of events with MadGraph 
at all. 
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Appendix B. Technical aspects of the MadGraph implementation 





rd>l dlllc Lei 


1V±gLU.vxl dpil 




Unit QtnVt MPPM 


Gauge 


\ 

ot 


aEWMl 


127.916 






Gf 


Gf 


0.000011664 


GeV A 




«s 


aS 


0.1184 








MZ 


91.1876 


GeV 


CKM matrix 




lamWS 


0.2253 







Aws 


AWS 


0.808 






Pws 


rhoWS 


0.132 






VWS 


etaWS 


0.341 




Fermion masses 


m c 


MC 


1.42 


GeV 




m t 


MT 


172.9 


GeV 




m b 


MB 


4.67 


GeV 




m e 


Me 


0.00051100 


GeV 






MM 


0.10566 


GeV 




m T 


MTA 


1.7768 


GeV 


Higgs masses 


m p > 


mrho 


125 


GeV 




m h < 


mhl 


200 


GeV 




m h " 


mh2 


150 


GeV 




m H ± 


mhc 


150 


GeV 



Table B.2.: The MCPM parameters and their names in the MadGraph implementation. 

Note that some parameters are zero in the MCPM (see last column), how- 
ever they can be set to other values in this model. This allows the correct 
simulation of SM background processes (see section 3.2). The CKM matrix 
is given in the Wolfenstein parametrisation \Y1\ section 7.8.1]. 
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Appendix C. 



Results of the validation 



C.l. Production cross-sections and decay widths 



The first check of the Monte-Carlo simulation consisted of the calculation of cross- 



sections and decay widths, see section 4.1. The results are listed in table C.l for the 



Higgs boson production and in table C.2 for the decay. 



Process 




Parameter set 


m H [GeV] 


CMadGraph [pb] 


otheo [pb] 


p + p ->■ 


h' + X 


1 


100 


13188.0 


±8.0 


13195.0 






2 


400 


75.382 ±0.048 


75.396 






3 


200 


1158.7 


±0.7 


1159.0 






4 


300 


244.44 


±0.15 


244.50 


p + p ->• 


ti' + X 


1 


100 


13199.0 


±8.0 


13195.0 






2 


200 


1159.0 


±0.7 


1159.0 






3 


150 


3284.7 


±2.0 


3284.5 






4 


250 


499.06 


±0.31 


499.10 


p + p — > 


H~ +X 


1 


150 


4454.7 


±2.6 


4454.5 






2 


250 


722.89 


±0.44 


722.87 






3 


200 


1644.4 


±1.0 


1644.4 






4 


300 


360.96 


±0.22 


360.95 



Table C.l.: The results of the systematic check of production cross-sections. Parameter 
sets and other details can be found in section 4.1.1| The deviations are very 
small, to say the least. 
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Appendix C. Results of the validation 



Process 


Parameter set 


m H [GeV] 


TMadGraph [GeV] 


Tthco [GeV] 


h! ->■ c + c 


1 


100 


5.8789 ±0.0072 


5.8789 




2 


400 


23.542 ±0.029 


23.542 




3 


200 


11.769 ±0.014 


11.769 




4 


300 


17.656 ±0.022 


17.656 


h" ->• c + c 


1 


100 


5.8837 ±0.0072 


5.8837 




2 


200 


11.771 ±0.014 


11.771 




3 


150 


8.8275 ±0.0107 


8.8275 




4 


250 


14.714 ±0.018 


14.714 


F~ ->• s + c 


1 


150 


8.8339 ±0.0108 


8.8339 




2 


250 


14.725 ±0.018 


14.725 




3 


200 


11.779 ±0.014 


11.779 




4 


300 


17.670 ±0.022 


17.670 


a) Hadronic decays. 


Process 


Parameter set 


m H [GeV] 


TMadGraph [keV] 


Ttheo [keV] 


fo' ->■ /X~ + /i + 


1 


100 


207.20 ± 0.25 


207.20 




2 


400 


828.82 ± 1.01 


828.82 




3 


200 


414.41 ± 0.51 


414.41 




4 


300 


621.61 ± 0.76 


621.61 




1 


100 


207.20 ± 0.25 


207.20 




2 


200 


414.41 ± 0.51 


414.41 




3 


150 


310.81 ± 0.38 


310.81 




4 


250 


518.01 ± 0.63 


518.01 


->■ /x~ + ^ 


1 


150 


310.81 ± 0.38 


310.81 


2 


250 


518.01 ± 0.63 


518.01 




3 


200 


414.41 ± 0.51 


414.41 




4 


300 


621.61 ± 0.76 


621.61 



(b) Leptonic Decays. 



Table C.2.: The results of the systematic check of total decay widths. Monte-Carlo 
results and theory agree perfectly. 
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C.2. Full processes: Additional plots 



C.2. Full processes: Additional plots 



Fig. C.l shows the distribution of pseudorapidity r\ in the lab frame for the n /x H 



channel. In fig. C.2 the resonance peak for the sc mode can be seen. In either case, the 



processes and parameters described in section 4.2.1 were used. 



0.006 
0.005 
0.004 
0.003 
0.002 
0.001 




Figure C.l.: Dependency of the /i n + cross-section on the pseudorapidity rj of the \i 
in the lab frame, yfs = 14 TeV. 
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Figure C.2.: Resonance in the sc channel with a Breit-Wigner fit. y/s = 14 TeV. 
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Appendix D. 

Results of the analysis of MCPM signatures 

D.l. Comparison of MCPM signal with background 



In section 5.1 the MCPM \i fi + signature was compared to the SM background. Fig. 



ID. II shows the cross-sections of MCPM and SM as a function of the invariant mass of 



the muons for ^/s = 7 TeV, both with and without cuts. In table D.l the size of the 



Higgs resonance peaks is compared to the background for a variety of settings. 
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D.l. Comparison of MCPM signal with background 



-v/s ITeVl 


Tnv Tnast; ran CP [OpVI 


r? • i ml>l 
"signal [P u \ 


u background LP u \ 




Ratio 

_L LCI LIU 


1 /I 
14 


1 /IO 1 
l^iZ - lOo 


Q OA Q 1 n - 2 

o.o4o • 1U 


0,01 


o.Uz l 


• 1U 




190 - 214 


3 334 • 1 0~ 2 


9 938 • 10 -1 


3.356 


• lO -2 




234 - 266 


1.444 . icr 2 


5.111 • lO -1 


2.826 


• lO" 2 




378 - 426 


2.149 • 10 -3 


8 894 • lO -2 


2.416 


• 10 -2 


7 


1 A9 1^8 
i^±z - loo 


9 4^8 . i n — 2 


1 998 


i .you 


. 1 n -2 




190 - 214 


8.086 • 10" 3 


4.078 • 10- 1 


1.983 


•lO- 2 




234 - 266 


3.159 • 10~ 3 


1.966 • lO" 1 


1.607 


• lO- 2 




378 - 426 


3.517- 10' 4 


3.045 • lO" 2 


1.155 


•lO" 2 


a) Raw data 


- no cuts applied. 










[TeV] 


Inv. mass range [GeV] 


^signal [pb] 


f background [pb] 




Ratio 


14 


142 - 158 


6.537-10- 2 


1.318 


4.959 


•lO- 2 




190 - 214 


2.783-10-2 


5.931 • 10- 1 


4.692 


•lO- 2 




234 - 266 


1.161 • 10-2 


2.665 • lO" 1 


4.356 


•lO" 2 




378 - 426 


1.589- 10- 3 


5.422 • lO" 2 


2.931 


•lO- 2 


7 


142 - 158 


2.010-10-2 


6.771 • 10" 1 


2.969 


• lO" 2 




190 - 214 


7.439 • 10- 3 


3.049 • 10- 1 


2.440 


•lO- 2 




234 - 266 


2.757- lO- 3 


1.301 • lO" 1 


2.119 


• lO" 2 




378 - 426 


3.218- 10" 4 


2.243 • lO" 2 


1.435 


•lO" 2 



(b) After applying cuts on rj and px- 



Table D.l.: The MCPM cross-sections at the Higgs resonance peaks. The Higgs masses 
are 150, 200, 250 and 400 GeV, respectively. 
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Appendix D. Results of the analysis of MCPM signatures 
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(a) Raw data - no cuts applied. 
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(b) After applying cuts on r] and pr- 

Figure D.I.: MCPM events compared to the SM background in the \x~ ' [i + channel as a 
function of the invariant mass of the muon pair. \fs = 7 TeV. 
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